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Abstract 

We discuss the prospects of observing leptonic CP-violation at a neutrino fac- 
tory in the context of the standard three neutrino scenario. If the large mixing 
angle MSW effect turned out to account for the solar neutrino deficit, we show 
that observing an asymmetry between the v e — > and v e — > oscillation 
probabilities would represent an exciting experimental challenge. We deter- 
mine the portion of the parameter space where an evidence could be found as 
a function of the intensity of the muon source and of the detector size for differ- 
ent baselines. We discuss the consequent requirements on a neutrino factory. 
We address the problems associated with asymmetries induced by the experi- 
mental apparatus and by matter effects. Finally, we emphasize the importance 
of measuring two CP-conjugated channels in order to precisely determine #13. 



1 Introduction 



An increasing evidence for the existence of neutrino oscillations has been gathered by 
several experiments. The standard theoretical interpretations of the data requires the 
neutrino flavour eigenstates to be superpositions of three or more mass eigenstates. This 
involves the possibility of a violation of CP associated with the physical phases of the 
mixing matrix. 

An evidence of CP-violation has been provided in the quark sector by different exper- 
iments, that have suggested a value of the CKM phase close to maximal. It would be 
extremely interesting to be able to investigate the CP-violation issue also in the lepton 
sector. Unfortunately, the upcoming long-baseline experiments will not help in this respect 
since they are not conceived for the purpose of comparing CP-conjugated transition. Note, 
however, that they could be affected by CP-violation in a relevant way. The oscillation 
probabilities they measure do in fact contain a CP- violating part. A detailed analysis [U 
has shown that the latter can constitute a relevant part of the full oscillation probability 
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in the v e *-* Un channel and could therefore affect the measurement of the mixing angle 

#13- 

An experimental study of CP-violation in neutrino oscillations would require a further 
generation of experiments offering the possibility of measuring oscillation probabilities in 
two CP-conjugated channels. Such possibility has been made concrete by the proposal 
of using the very intense muon sources that are currently being designed as part of the 
muon collider projects to produce an high-intensity neutrino beam j^, ^, |]]. As a pure 
source of both neutrinos and antineutrinos with well known initial flux, such a "neutrino 
factory" would be an ideal framework for studying leptonic CP-violation. The possibility 
of exploiting a neutrino factory for a measurement of CP-violation has been first considered 
in 0, where the capabilities of a reference set-up have been studied in different points 
of the neutrino parameter space. An update of that analysis can be found in |J]. The 
purpose of this paper is to explore the full parameter space and to show in which part of 
it the detection of a CP-violating effect would be possible. The portion of the parameter 
space that could be covered will be determined as a function of the intensity of the muon 
source and of the size and efficiency of the detector for two different baselines. This 
will make clear which values of the factory parameters are necessary in order to cover 
a given portion of the parameter space. We address these issues in the context of three 
neutrino oscillations. Analysis of CP-violation in four neutrino oscillations can be found 

in mm. 

The paper is organized as follows. In Section ^ we introduce the usual notations for 
the neutrino masses and mixings and we describe the experimental set-up. In Section |3] we 
discuss the problems involved with a realistic measurement of CP-violation. In particular 
we describe how we cope with the asymmetry induced by the experimental apparatus and 
by matter effects. The systematic and statistical errors associated with the measurement 
are also discussed. In Section |I| we determine the optimal values of some parameters of 
the experiment and we present the results about its capabilities for different experimental 
configurations. In Section |5] we discuss the results and we conclude. 

2 Framework 

2.1 Masses, mixings and probabilities 

In this paper we consider the possibility of measuring a violation of CP in the leptonic 
sector by means of a neutrino factory Q. We address this issue in the context of three 
neutrino oscillations. The leptonic charged current is then given by 

err'Ve, = ea^UijVj, (1) 

where Uij is the 3x3 unitary lepton mixing matrix and ej, i/j, i = 1 . . . 3 are the left-handed 
charged lepton and neutrino mass eigenstates. By convention we order the neutrino mass 
eigenstates in such a way that m Ul < m V2 and Am^ < | Am^j | , | Am| 2 1 , where we denote 
as usual Am| = m^. — . The common sign of Am^ and Am| 2 is not determined at 
present. 



2 



We use the standard parameterization for U, 



(C12C13 S12C13 si 3 e \ 

-S12C23 - ci2S23Si 3 e J<5 C12C23 - si 2 s 2 3Sne t5 S23C13 , (2) 
S12S23 - Ci 2 C23Sl3e J<5 -C12S23 - Si2C23Sne lS C23C13 / 

where possible Majorana phases are omitted because they do not enter the oscillation 
formulae in the safe approximation in which lepton number violating oscillations are ne- 
glectedQ. In that approximation, the v ei — > u ei , u ei — > v ej oscillation probabilities in 



vacuum are 

P(v ei -> v e . ) = P C p {v ei ~^v ej )+ Pop (v ei -> v £] ) (3a) 

Pipe, Ve,) = PCP iy e . -> V e .) - Pqp (u e . U ej ), (3b) 

where 

Pcp (y ei -^v e .) = Sij - 4 Re J( l 2 sin 2 Ai 2 - 4 Re J d 2 \ sin 2 A 23 - 4 Re sin 2 A 3i , (4a) 

P&p(y ei f G .) = -8ciijJsinAi2sinA23sinA 3 i, (4b) 

J kh = u e t u k Ut k e 3 U ejVh Ut hei , A kh = Am 2 kh L/4E, = Y.k e ijk and 

8 J = cos 0i3 sin(26>i 3 ) sin(26»i 2 ) sin(26>2 3 ) sin 5. (5) 



The ranges we use for the neutrino masses and mixings are those obtained by the 
standard fits of the atmospheric and solar neutrino data H |9[ [l0|, [0]] taking into account 
the constraints given by the CHOOZ experiment |j~2[| . Among the three possible standard 
solutions of the solar neutrino problem we only consider in detail the "large angle" one. 
In the "small angle" solution, as well as in the "vacuum" solution, CP-violation effects are 
in fact too small to be detectable, since they are suppressed either by the small value of 
012 (small angle solution) or by the small value of Am 2 2 (vacuum solution). 

2.2 Experimental set-up 

The possibility of using the straight section of a high intensity muon storage ring as a 
neutrino factory has been emphasized in ||]. Being a pure and high intensity source of 
both neutrinos and antineutrinos with well known initial flux, a neutrino factory would give 
the possibility of comparing the oscillation probabilities in two CP-conjugated channels: 
v ei — > v e .jv ei — > v ei (i ^ j). Therefore, it would be the ideal framework for studying 
leptonic CP-violation. Among the possible channels made available by a neutrino factory, 
the most suitable for studying CP-violation are v e — > u^/Uf, — > and the T-conjugated 
— > u e /i>^ — > v e . In fact, in these channels the CP-violating part of the oscillation 
probability is not submerged by the CP-conserving part (as it is for the <-> v T channels), 
so that large asymmetries between the CP-conjugated channels can arise [y. We consider 
here in particular the v e — > v^/v e — > channels, since telling /j> + from in a large high- 
density detector is easier than telling e + from ||. 

1 We prefer to think of three Majorana neutrinos but in the mentioned approximation what follows 
apply to the case of three purely Dirac neutrinos as well. 
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The v e — > v„ [p e — > v^) oscillation probability can be measured as follows: the electron 
neutrinos (antineutrinos) are produced by the decay of iV„+ positive muons (N^- negative 
muons) in the straight section of the storage ring pointing to the detector; then the Va 
(y^) produced by the oscillation can be detected by their charged current interactions in 
the detector. The number n(z^) (n(i^)) of "observed" (pn) is then given by 

n(u,) = N^ + N kT ^i-% / /(£ )jP ™ (£) (6a) 

in9AJ". E^ /'Bmax _ J TP 

n{9») = %-N kT —^^ / f(E)P^(E) — , (6b) 

f_l J -Emin M 

where 



f(E) = g^E/E^a^/E^-iE), (7a) 
f(E) = g Pe (E/E^)(a^(E)/E^ + (E) (7b) 

are weight functions taking into account the energy spectrum (normalized to 1) of the elec- 
tron neutrinos (antineutrinos) produced in the fi + decay, g Ue ^(E/ E^), the charged 
current cross section per nucleon, a u ^ p ^{E), and the efficiency for the detection of fj,~ 
e^-( M +)(-EO0- P^-^u (E) (Pt%->v (E)) 1S tne oscillation probability for neutrinos 
travelling inside the earth taking into account matter effects. Finally A^t is the si ze of 
the detector in kilotons and JV„± is the number of "useful" fj, decays, namely the number 
of decays occurring in the straight section of the storage ring pointing to the detector. In 
the numerical calculations we use 

3, e (z)= 5 p e (x) = 12x 2 (l-x), (8a) 
a Vft (E) = 0.67 • 1(T 38 £ cm 2 /GeV, {E) = 0.34 • 10~ 38 £ cm 2 /GeV, (8b) 

and e^± (E) = 30% for E > 5 GeV, (8c) 

so that f(E)/f(E) = 2 (independently of the energy). We also apply a lower cut -Emm = 
5 GeV on the neutrino energies in order to have a good detection efficiency in all the energy 
rangef^. 



3 Measuring CP-violation 

3.1 Experimental apparatus and matter asymmetries 

The asymmetry between the number of i>a and 9^ events seen in the detector, normalized 
to the initial number of muon decays is given by 

atot n{v,)/N^ + + n(* M )/JV ' 1 ' 



2 We neglect here the finite resolution of the detector. 

3 The numerical results we obtain agree with those obtained in [g] if the efficiency is set to 100% and 
no lower cut is used. 
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An "intrinsic" leptonic CP-violation associated with a non-vanishing phase 5 in (||) would 
contribute to give a non- vanishing atot- However, even if there were not any intrinsic CP- 
violation, citot would still be non-vanishing due to the CP-asymmetry of the experimental 
apparatus, / ^= f in (g), and due to matter effects, that affect neutrinos and antineutrinos 
in different ways. In other words, atot is an estimator of 

J P^(E)f(E) dE-J P^(E)f(E) dE 
atot JP^(E)f(E)dE + JP^(E)f(E)dE J 1 j 

that can be non- vanishing even when 5 = due to / ^ f and P™^ u (E, L) ^ Ppl^v (E, L). 

Getting rid of the experimental apparatus asymmetry is easy. It is sufficient to weight 
e.g. each v event with f(E)/f(E), where E is the energy of the event, namely to replace 
n(^ At ) in (^) with 

n '(p ) = y ik^l, (11) 
where Ei is the energy of the i-th event. In other words, instead of atot we consider 

that estimates a quantity independent of the experimental CP-violation, 

_ JXgWg - P^(E))f(E) dE 

aCP+m !{ prn^ {E)+ p^ {E))f{E)dE - I ) 

We remark, however, that the subtraction of the experimental asymmetry introduces an 
unavoidable systematic error due to the uncertainties in the knowledge of the flux, cross 
section and efficiency ratios in f(E)/f(E). 

In our numerical calculations n'(pn) is simply given by 

n'(^)=2n(^), (14) 



since f(E)/f(E) = 2 independently of the neutrino energy E. In this case, using (12) 
is equivalent to using a different normalization for n(i^) and n(u^) in (0). However, if 



f(E)/f(E) had a small dependence on the energy due e.g. to e ll -(E)/e lx +(E), using (12) 
would allow a more accurate subtraction of the experimental asymmetry. 



The quantity dcp+ m in (12) still contains a pollution due to the matter asymmetry. 
In principle the best attitude towards such a pollution would be trying to experimentally 
disentangle matter from genuine CP-violation effects. A detailed study of the feasibility 
of this possibility is beyond the aim of this paper. Here we would rather cope with matter 
effects by calculating them [||. An experimental confirmation of the reliability of such a 
theoretical calculation, besides being important in itself, would be welcome in order to 
support this method. Other theoretical approaches to the issue of matter effects can be 



found in 113 
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First of all, we define the matter asymmetry in absence of CP-violation as 



f(P^(E) - P™j=»(E))f{E)dE 
f(P™!l=°(E)+P™if(E))f(E)dE : 



III.:): (»,,-. . r> ,ll.:). . . , • 

where P™'!^,®, P^'J^® are the oscillation probabilities calculated for 5 = but taking into 
account matter effects. Then, we define the CP-asymmetry in absence of matter effects as 

= JiPj^jE) - P Pe ^(E))f(E) dE = J Pqp (E)f(E) dE 
flCP J(P Ue ^(E)+P^(E))f(E)dE jP CP (E)f(E)dE , [ ' 

where P Ue ^ Uti , P Pe ^p are the oscillation probabilities computed in vacuum but taking 
into account CP-violation effects. Now, despite acp+ m is not simply given by the sum of 
acp and a m , it turns out that for practical purposes the relation 

OCP+m - «CP + «m (17) 

is a very good approximation if a m is not too largeQ. In any case, the error one makes 
recovering acp through (|l7|), acp = acp+ m -«mi is smaller than the uncertainties on a m 
itself. The latter therefore constitute the main source of systematic error from matter 
effects. In conclusion, as a measure of CP-violation we use the quantity 

_ n(^)/A>-n'(^)/iy 

Being dcp an estimator of acp in (|i~6|), measuring a non- vanishing value of dcp would 
indicate a genuine leptonic CP-violation. 



3.2 Systematic and statistical errors 

In this Section we discuss the systematic and statistical errors involved in a measurement 
of acp through acp- There are different sources of systematic errors in acp- One is given 
by the uncertainties on the experimental asymmetries, g Ue /gp e , ^val a Oa an d ^-/e^+j that 
enter acp through ///in (|TT|). Another one is given by the uncertainties on the matter 
asymmetry, in turn due to the uncertainties on the neutrino mass and mixing parameters 
we use to calculate it. Here we concentrate on the latter, assuming that the former will 
be made small enough by experimental determinations of ///. 

Whatever is the strategy used towards the matter asymmetry, it is clear that the 
smaller the matter effects are the cleaner a possible measurement of a genuine violation 
of CP would be. Since the matter asymmetry grows with the baseline faster than the 
CP-asymmetry, a m constitutes the main limitation to the length of the baseline of an 
experiment aiming at measuring a violation of CP. In order to illustrate and make quan- 
titative the latter point we can use the approximation one gets for , P^^ p ^ in the 

4 This is essentially because the corrections to (^) arise only at second order in the asymmetries and is 
confirmed by a numerical analysis flf. 
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limit Am 2 ^ = 



= Bin 2 23 sin 2 2 g 13 ^ 2 ( A 32cf) (19a) 



^r% ; sin^^ Sin2(A32C - 2) , (101)) 



where 



c± = (lip x ) 2 ±4x sin 2 6»i3 ~ (1=fx) 2 , x = ^— g- (20) 



:!2 



and V is the matter induced potential. We then get 

; + 7 )/c+ - sm^(A 32 c_ / 
" sin 2 (A 32 cf)/c + + S in 2 (A 32 cL /2 )/c_ ~ 12 E * 

for the matter asymmetry before integration on the energy range, or 



sin 2 (A 32 4H/c + - sin 2 (A 32 cr)/c_ ^ VAm 2 32 L 2 , > ; 2 



2sin z # 13 ) (21) 



a L 2 (km) Ami , . 

°^°- 7 - 10 g (GeV)3.10-^ - <22> 

The expansion in (21) holds when L < tt/V and A 32 < ir/2. When these conditions 
are not fulfilled the matter asymmetry is large or quickly oscillating and therefore out 



of control. However, we stress that the main corrections to (22) come from the approx- 
imation (|19|), particularly rough when #i 3 is very small, rather than from the expansion 
in (^). Nonetheless, Eq. (|22^ ) is sufficient for our illustrative purposes (exact formulae are 
used in all numerical calculations). When compared with the behavior of acp with L and 
E, 

a C p oc sin A 2 i oc ^, (23) 

Eq. (p2^) shows that |a m /acp| grows with the baseline length L (as far as L < tt/V and 
A32 < tt/2). Very long baselines are therefore more suitable for the study of matter effects 
than for CP-violation. We will come back to this point in Section ^J. 

Before describing into detail how we work out the systematic error due to matter effects, 
we remark that Eq. ( p2] ) , besides illustrating the linear dependence of a m on Am 2 2 , makes 



also explicit its dependence on the sign of Am 2 2 . The knowledge of that sign is therefore 
essential in order to properly subtract matter effects and will be assumed here. A detailed 
analysis of how that sign could be determined through a measurement of matter effects is 
beyond the scope of this paper. 

We estimate the systematic error <5 sys tacp on a CP due to the uncertainty on a m by 
scanning the ranges of the parameters on which a m depends. 5 sy stacp is then obtained as 
the half-width of the range of values assumed by ti m . As for the Am^, 6q 2 , # 23 ranges, 
less important in the determination of <5 sys tacp, we assume, to be conservative, that they 
will not differ very much from the present ones ^ pi |ll[ : 

Amli = (0.2 - 3.0) • 10~ 4 eV 2 , sin 2 29 12 = (0.6 - 0.96), sin 2 26 23 > 0.8 (24) 
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(all at 99% CL; the Amfi, #12 values correspond of course to the large mixing angle 
solution). It is less clear which ranges should be used for the more crucial Am^ and #13. 
The next generation of long-baseline experiments will be able to provide information on the 
latter parameters and will reduce the uncertainties on them. However, it has been shown 
in that, for favourable values of Am^ in the large angle range in (|24|) and for values of 
Am 2 2 and #13 within the sensitivity of an experiment like e.g. MINOS, the CP-violating 
part of the oscillation probability could constitute up to 30-40% of the total oscillation 
probability. This in turn means that, unless some compelling evidence excluding the large 
angle solution is found, a systematic uncertainty up to 30-40% could affect a sin 2 26*13 
determination made measuring a single oscillation probability. A precise determination 
of 0i3 would on the contrary require an experiment able to measure both the u e *-* Va 
and 9 e *-* oscillation probabilities, namely able to measure CP-violation. A neutrino 
factory would provide this possibility. Once both P Ve -+ u and P Pe _» 5 were measured, the 
CP-conserving part of the probability could be recovered as Pqp = {P v& -± v + Pp e _,.p )/2. 
This would not completely remove the ambiguities associated to the phase S, since Pqp 
depends on 5 as well. An estimate of 5 through acp would therefore be also welcome in 
order to get rid of such ambiguities. Here, for the purpose of estimating <5 sy stacp we will 
assume that a determination of Am| 2 and sin 2 2#i3 with a precision higher than 20% will 
become available. 

The discussion of statistical uncertainties is straightforward. In order to exclude the 
possibility that a measurement of a non-vanishing acp is due to a statistical fluctuation, 
the measured value must be larger than n CT <!> s tatacP; where <5 s tatacp is the "lcr" statistical 
error on acp in absence of CP-violation and n a is the number of standard deviations we 
require. Since in absence of CP-violation the expectations of n(v^) / and n'(^)/iV„- 
are equal, we get 

1 1X1/2 



<WacP = 77-7 — rr + . , rr , (25) 
V4(n(z^)) 4(n(i/ M ))/ 

where (n(z/„)) and {n{i>^)) are the expected number of Vn and interactions seen in the 
detector. When deciding whether it would be possible to measure a violation of CP in a 
given point of parameter space, we will assume as usual that the measured value of acp 
coincides with the value of acp expected in that point. 

Finally, let us shortly discuss the contribution of the background to the statistical 
error. According to the present estimates, the main source of background is due to charm 



production in the charged and neutral current neutrino interactions in the detector [14|. 
For example, the charged current background is due to: 

_ noose. CC int. / _ N c-decay , i. , . 

M -> Vfl * "fi ► (/X ) lost + C > (M ) found- (26) 

This background only affects the signal corresponding to the v e — ► oscillation: 

_ _ osc. _ CC int. 4. tnn\ 

H -> v e ► Ufj, > /X . (27) 

The background subtraction corresponds to a replacement 

n(z^) -> niufj) - (n(P M ))back (28) 
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in Eq. (|18|), Such a subtraction introduces a further source of statistical error. Using the 
estimate jL4| 



(n(^))back ~ 10-V(^) P=1 , (29) 

where n{v^) p= ^ stands for the number of 9^ interactions that would be seen if all the 
initial 9 e oscillated into 9^, we get the complete expression for the statistical error we will 
use: 



J 1 10- 5 (n(9, Jp= li 

4<n(^)> 4(n(^)) 2 



Note, however, that in a long-baseline experiment with sin 2 A32 = O (1) the background 
contribution in (|30| ) is not negligible only for sin 2 2^13 < 10 -5 . 

4 Capability of the neutrino factory 

In this Section we show in which part of the parameter space an evidence of leptonic 
CP-violation could be obtained as a function of the relevant experimental parameters. In 
order to obtain an evidence for a non-vanishing 

f(P Ue ^(E)- P Pe ^(E))f(E)dE fPop(E)f(E)dE 



a ° P I(P,[Jl(E)+P p [J^E))f(E)dE jP CP {E)f{E)dE' ® 
a value of acp larger than the systematic and statistical errors on it should be measured: 

|6cp| > 4 y stacp, |acp| > ^Atatacp- (31) 



In the previous equation, <5 sys tacp is calculated as described in Section p,2| , <5 s tatacp is the 
la statistical error given by ( pOj ) and n a is the number of standard deviations corresponding 
to the CL we want to achieve (we will use n a = 5 as a reference value). 
The experimental parameters we consider are 

• the size of the detector in kT, N\^; 

• the total number of useful muon decays = N„+ + iV„- ; 

• the ratio of \x r versus \i~ decays iV„+ / N„- for a given total number of muon decays 

• the energy of the stored muons En and the lower and upper cuts on the neutrino 
energies, E min and E max ; 

• the baseline L. 

First, we briefly identify a convenient choice of N^+fN^- and E^, E m \ a , E max . Then, 
using that choice, we show which is the portion of parameter space where the asymmetry 
could be measurable as a function of the product NuN^ and for two choices of the baseline, 
L = 732 km and L = 3000 km. 
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The dependence of the statistical error on and N^- is simple. Neglecting the 
small background contribution and using Eqs. (|8|) we find 

<W«CP oc (— !— + Tf— — ) , ( 32 ) 

where ov M /o"p M ~ 2. When A r At + = JV„- the ^ events are less than the events because of 
the smaller i/„ cross section. Therefore, their contribution to the statistical error is larger. 
For a given total number of available muon decays, = JV„+ + A r At - , the statistical error 
can be made smaller by choosing a ratio N^+fN^- smaller than one. Here we therefore 
use the / N^- ratio obtained by minimizing ( |32[ ) , namely 

A> = ' iV M ^ 0.4 (33a) 




iV = ^-l* ^ - 0.6 N„. (33b) 



Unfortunately, using ( |33| ) instead of A^ M + = AL- only improves <5stat«CP by a few percents. 

Let us now discuss the energy parameters. First of all, we use E m \ n = 5 GeV as a 
lower cut on the neutrino energies in order to have an efficient muon detection in all the 
energy range. An upper cut is also necessary since acp/5 s tatOcP oc 1/V E m3iX when E max 
is sufficiently high |5). The upper cut will of course be lower than the muon energy in 
the storage ring. The low energy neutrino rates do not strongly depend on the value 
of Efj, > E max , as far as E^ is in the range presently discussed 10 GeV < E^ < 50 GeV. 
However, this dependence is not negligible either. For example, the rate of charged current 
interactions of low energy neutrinos in the detector is inversely proportional to E^, as it 
can easily seen by using Eq. (|8|). 

The optimization of E max and E^ is shown in Fig. |l] for the two baselines we consider, 
L = 732 km (Fig. |l|a) and L = 3000 km (Fig. |l]b). The values of E max and that 
maximize |acp/<5 s tatacp| are shown by the black spots in the E^-E^j 'E m i n planes of Fig. [l|. 
The three contour lines around the spots surround the regions where |<2cp/^stat^cp| is 
larger than 80%, 60%, 40% of the optimal value respectively. The figures show that there 
is a certain freedom in choosing the values of E max and E^. In what follows we will use the 
optimal values E max = E^ = 20 GeV in the L = 732 km case and E max = E^ = 40 GeV 
in the L = 3000 km one. However, using e.g. E max = 20 GeV and E^ = 50 GeV for 
L = 732 km would involve a reduction of |acp/^statOcp| of 20% only. The dependence of 
these results on the neutrino parameters is negligible. 



Having described the values of N^+jN^- E^, E m i n , E max we will use, we now move 
to the discussion of what amount of parameter space could be covered by using a total 
number of A^ = A^+ + N^- useful muon decays and a kT detector. First of all we 
describe the parameter space. The neutrino parameters are #12, O23, #13, Am^, Am| 2 , S. 
The constraints on them we use are those obtained by the standard fits of the atmospheric 
and solar neutrino data || [| |To|, |Ti|] taking into account the constraints given by the 
CHOOZ experiment [Oj. Since the uncertanties on #12 and #23 do not significantly affect 
the results, we set sin 2 2#i2 = 0.8 || and sin 2 26*23 = 1- The dependence of the results on 
5 is trivial since acp is to a good approximation linear in sin<5. This is because Pqp is 
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(a): L = 730 Km (b): L = 3000 Km 




E max (GeV) E max (GeV) 

Figure 1: Optimization of E max and E^/E max in the case of a long baseline of L = 732 km, 
(a), and of a very long baseline of L = 3000 km, (b). The black spots indicate the values 
ofE max and E^j E max that give the best relative statistical error. The three contour lines 
surround the regions where the relative statystical error is 80%, 60%, 40% of the optimal 
value respectively. 

strictly linear in sin 5 and Pcp depends on 5 only through sub-leading Am^i terms[]. We 
are therefore led to consider a three-dimensional parameter space described by 

10~ 3 eV 2 < |Am| 2 | < 10~ 2 eV 2 § | (34a) 
2 • 10~ 5 eV 2 < Am^ < 3 • 10" 4 eV 2 @ (34b) 
013 < CHOOZ and 3-i/ fit limits 0. (34c) 

Since the CP-violating effects crucially depend on Am^, we plot our results in the 
sin 2 2#i3-Am 2 . 1 plane for different values of |Am 2 2 |. We first show and discuss the re- 
sults for the L = 732 km baseline. 

4.1 The L = 732 km baseline 

In Fig. H we show the portion of the sin 2 2#i3-Am 21 parameter space where an evidence 
of CP-violation could be obtained for three different values of | Am 2 2 1 : the present central 
value, |Am| 2 | = 3 • 10 -3 eV 2 (Fig. |2|b in the center), and two values at the borders of the 
presently allowed region, 

|Am| 2 | = 10~ 3 eV 2 (Fig. |a) and |Am§ 2 | = 10" 2 eV 2 (Fig. |c). 
The rectangular windows inside the light shadowed areas represent the part of parameter 
space allowed by the constraints (|4|). The dark shadowed parts of the plots represent 
the portion of parameter space where an evidence of CP-violation could not be obtained 
because the systematic error exceeds the asymmetry itself even for maximal CP-violation. 
We see that the systematic error due to matter effects does not represent a serious problem 
for the L = 732 km baseline, especially when |Am 2 2 | is low. Note that the size of this 

5 Moreover, the CP-conserving quantity Pcp depends on S only through cos S, so that the dependence 
on S is particularly mild when the CP-violation is large. 
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Figure 2: Capability of a neutrino factory using a baseline of L = 732 km in the sin 2 2^13- 
Am^ plane for |Am§ 2 | = 10~ 3 eV 2 , (a), |Am 2 2 | = 3- 10~ 3 eV 2 ; (b), and |Am| 2 | = 
10 _2 eV 2 ; (c). The rectangular windows inside the light- shadowed areas correspond to the 
constraints $32)- In the dark-shadowed regions the systematic error on acp due to matter 
effects exceeds \acp\- The size of these regions depends on the precision of the future 
determination of |Am| 2 | and #13. In order to cover the portion of parameter space above 
the solid line corresponding to a given value of V , values of N^, N^t an d e^± such that 
N^N kT /(n a /5) 2 (e^± /30%) sin 2 S > V should be used. 



systematic error depends on the precision of the future determinations of |Am 2 2 | and #13. 
If precisions higher than those assumed here could be achieved, the dark-shadowed regions 
in Fig. ^ would be smaller. Fig. ^ shows that the relative systematic error |<5 sys ttJcp/acp| 
grows with |Am 2 2 |. This is because acp is almost independent of Am| 2 when #13 is not 
too small whereas |<5 syst acp|, roughly proportional to |a m |, grows with |Am| 2 |, as 
illustrated by Eq. (g|). 

The region of the parameter space where the statistics is sufficient to exclude a fluc- 
tuation with a CL corresponding to n a standard deviations is shown in the plots of Fig. [2| 
for different values of the number of useful muon decays Nn, of the size of the detector 
in kT, iV^T, and of the size of CP-violation, | sin5\ (regions above the solid lines). Since 
acp — acp(S = vr/2) • sin 5 and since <5 s t a tOcp oc l/y/N^N^, the regions shown actually 
depend on the combination 

v = ^JWso%- sm 6 (35) 

only. For a 30% efficiency, maximal CP-violation and a CL corresponding to the reference 
value n a = 5, V is simply the product of and A^x- The three solid lines in Figs. |^ 
(one in Fig. ^a) correspond to three possible values of V: the smallest one, V = 4 • 10 21 , 
corresponds (for sin 2 5 = 1 and n a = 5) to the reference set-up considered in ||, a 10 
kT detector and 2 • 10 20 useful muon decays in both the CP-conjugated channels. The 
intermediate value, V = 4 ■ 10 22 , corresponds e.g. to an improvement of one order of 
magnitude in the intensity of the muon source. The highest value, V = 4 • 10 23 , requires 
a further improvement, e.g. a 30 kT detector with very high efficiency. Some of these 
options are currently under discussion [15]. 
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Figure 3: Contour lines of\acp(b~ = 7r/2) | ~ |acp/sin5| in the sin 2 26 ) i3-Am 21 plane for 
\Aml 2 \ = HT 3 eV 2 , foj, |Am| 2 | = 3- l(T 3 eV 2 , (b), and |Am| 2 | = 10 _2 eV 2 , (cj, in the 
case of a L = 732 fcm baseline. In (c) the contour lines inside the darker region correspond 
to a 70% asymmetry. 



Fig. |2] clearly shows that the sensitivity lines strongly depend on the value of | Am 2 2 1 . 
This is because the statistics, proportional to sin 2 A32, grows with |Am 2 2 |, whereas the 
asymmetry is almost independent of Am| 2 for #13 not too small Q. A |Am| 2 | value 
around the present central value would allow to cover a major part of the parameter space 
with the highest value of V . If |Am| 2 | were at the top of the presently allowed range, 
that value would allow to cover essentially all the parameter space, while the intermediate 
value of V in this case would be enough to cover a good half of it. At the bottom of 
the |Am 2 2 | range getting an evidence of leptonic CP-violation would only be possible in 
a small part of the parameter space even with the highest neutrino fluxes, unless a longer 
baseline is used. This possibility will be considered in the next Subsection. 

Values of Am 21 above the range (|34|b) are allowed in "non-standard" analysis of solar 
neutrino data where the possibility of an unknown systematic uncertainty in one of the 
solar experiments is taken into consideration fl6|| . Fig. ^ shows that relatively low values 
of V could cover the corresponding region. 

For completeness, we show in Fig. ||three contour plots of |acp(<5 = vr/2) | ~ |acp/ sin5| 
in the same sections of the three-dimensional parameter space as in Fig. ^. The asymmetry 
values corresponding to the contour lines are specified in the legends. The structure of 
the plots is easily explained: as far as the Am 21 effects in the CP-conserving part of the 
oscillation probability can be neglected, 

Pqp oc sin 26*13 and Pep oc sin 2 26*13, (36) 

so that the CP-asymmetry gets larger when #13 gets smaller, acp oc 1/ sin 2#i3 Q. How- 
ever, at some point the smallness of #13 makes the Am^ effects in Pqp important so that 
acp can vanish when #13 — > as it should. 

Large asymmetries are therefore associated with small #13 and with small oscillation 
rates. One can then wonder whether having a large asymmetry would not represent 
a disadvantage from the point of view of statistics. Nevertheless, this is not the case. 
From (|3g|) we see in fact that, as far as we can set Am 21 = in Pep, acp/^statocp 
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Figure 4: Same as in Fig. |] for a baseline of L = 3000 km. 



Pqp I V Pcp is actually independent of #13. This explains why the sensitivity lines in 
Figs. [2] are approximately horizontal as far as #13 is not too small. 



4.2 Very long baselines 

We have seen in the previous Subsection that the possibility of covering a large region of 
the parameter space with a L = 732 km baseline crucially depends on |Am| 2 |- A better 
determination of |Am| 2 | will come in the next years from the upcoming long-baseline 
experiments. If |Am| 2 | turned out to be close or above the present central value and if 
the fluxes under discussion could be achieved, the L = 732 km baseline could be enough. 
If, on the contrary, |Am 2 2 | turned out to be low or the CP-violation size |sin 2 5| were 
small a longer baseline would be necessary. The reason why a longer baseline helps from 
the point of view of overcoming the statistical error is that the asymmetry grows with L, 
whereas the statistics, and therefore <5 s tatflCP! is approximately independent of L. This 
holds however only as far as |A 32 | <C vr/2, or L < 400 km^ GeV)(3 • 10 -3 eV 2 /| Am| 2 |). 

In this Subsection we consider the possibility of a L = 3000 km baseline. One can won- 
der whether having an even longer baseline could not to be better. However, this is not nec- 
essarily the case. First of all, sin 2 A3 2 reaches its maximum when L ~ 400 km E v { GeV)(3- 
10 -3 eV 2 /| Am 2 2 |), so that for longer baselines the statistics decreases significantly. More- 
over, matter effects grow with L faster than the asymmetry, as pointed out in Section |3.2| . 
This is also confirmed by a comparison of Fig. ^ with Fig. |], that only differs from Fig. [2| 
for the length of the baseline, L = 3000 km instead of L = 732 km. We see that the 
dark-shadowed regions where the matter effect uncertainties are too large are larger in the 
present L = 3000 km case. Therefore, for longer baselines matter effects would represent 
a major problem. Note, however, that these uncertainties depend on the precision of the 
future determinations of |Am3 2 | and #13 and they could therefore be smaller than what 



estimated here (see Section 3.2 



We see from Fig. [| that using a L = 3000 km baseline and the highest value of V 
would allow to cover about 2/3 of the Am 21 parameter space (for sin 2 26>i 3 > 2 • 10~ 3 ) in 
the most unfavourable case |Am 2 2 | = 10 -3 eV 2 and essentially all of it for larger values 
of |Am| 2 |. For such values, the intermediate value of V would be enough to cover a good 
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Figure 5: Same as in Fig. [3|/or a baseline of L = 3000 km. 



half of the Am|i parameter space. 

Finally Fig. || shows the contour plot of \acp(S = tt/2)\ ~ |acp/sin5| for the L = 
3000 km baseline. 



5 Discussion and conclusions 

The outcomes of the upcoming neutrino experiments will play a crucial role in assessing 
whether it could be worth making an effort towards a CP-violation measurement. First 
of all, they will help to understand whether a light sterile neutrino exists and mixes 
significantly with the active ones. If this were the case, the detection of leptonic CP- 
violation would be within the capabilities of an intermediate-baseline experiment jl], [|. 
If, on the contrary, the standard three neutrino scenario were confirmed, the size and the 
possibility of measuring CP-violation would crucially depend on the mechanism accounting 
for the solar neutrino deficit. If the large mixing angle solution were ruled out, there would 
be no hope of measuring leptonic CP-violation through neutrino oscillations. In the case 
of the small angle solution, this can be seen using Figs. ^ and |j. In fact, those figures also 
apply to the small angle solution provided that the parameter V is interpreted as 

r= MWm%- sm 6sm 2 * 12 ' (37) 

where 6\2 is now constrained to be within the small angle range. Therefore, in order to get 
same coverage in the sin 2 2#i3-Am 2 1 plane as in the large angle case, it would be necessary 
to increase • by a factor 1/ sin 2 2#i2 ~ 200]^| If, on the other hand, the large mixing 
angle solution were preferred by the data, the measurement of a leptonic CP-violation in 
the v e <-> valise «-> Vn channel would represent an exciting experimental challenge. 

In the latter case, a precise assessment of the capabilities of a neutrino factory, as well as 
the determination of the best experimental configuration (essentially the baseline length) 
would strongly depend on the value of two parameters: the intensity of the muon source 



6 Moreover, that would not still be enough since the Am^ range of the small angle solution is approxi- 
mately one order of magnitude lower than the Am^ range of the large angle solution. 
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and the precise value of |Am| 2 | which will be provided by the long-baseline experiments. 
Figs. [2] and || show what portion of parameter space could be covered for different values 
of these two parameters and for two possible baselines, L = 732 km (Fig. |2|) and L = 
3000 km (Fig. Q). If |Am 2 2 | turned out to be in the upper part of its present range, 
|Am| 2 | ~ 10~ 2 eV 2 , and very high-intensity sources could be achieved, the L = 732km 
baseline could be enough to cover a relevant part of the parameter space. On the other 
hand, if |Am| 2 | turned out to be lower, using an even longer baseline, L ~ 3000 km, could 
be necessary. Using baselines longer than 3000 km would require knowing matter effects 
with a higher precision than that assumed in this paper. Note that, since the relative 
systematic error due to matter effects grows both with L and |Am 2 2 | (as opposed to the 
statistical one, that decreases both with L and lAm^l) the very long baseline option could 
turn out to be less appropriate for the higher values of lAm^l- 

Finally, a neutrino factory would give the possibility of measuring #13 with a high 
accuracy. This is not only due to the pure and high-intensity neutrino flux, but also due 
to the possibility of measuring both the v e «-> i/„ and v e «-> channels. A determination 
of sin 2 2#i3 made by measuring a single oscillation probability only could in fact be affected 
by a systematic uncertainty up to 30% associated with CP-violation effects (]J (see also 
Figs. H and ||). On the contrary, measuring both P Va ^ v and Pp e ^p M would allow to get rid 
of the CP-violating part of the probability Pqp by computing (P Ve ^ v +Pp B -+v )/2 = Pqp 
and to get rid of the further dependence of Pqp on 5 by computing the CP-asymmetry 
ocp- 
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